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Abstract 

In this paper, the synchronization in a hyper-network of coupled dynamical systems is investigated for the 
first time. An evolving hyper-network model is proposed for better describing some complex systems. A 
concept of joint degree is introduced, and the evolving mechanism of hyper-network is given with respect to 
the joint degree. The hyper-degree distribution of the proposed evolving hyper-network is derived based on 
a rate equation method and obeys a power law distribution. Furthermore, the synchronization in a hyper- 
network of coupled dynamical systems is investigated for the first time. By calculating the joint degree 
matrix, several simple yet useful synchronization criteria are obtained and illustrated by several numerical 
examples. 
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1. Introduction 

Complex networks P, 0, 0, HI HI have been used to model the real complex systems with a large number of 
interacting individuals, such as Internet Q, world trade web 0, metabolic networks Q, coauthor networks 
9], and so on. The individuals in complex systems are denoted by nodes in complex networks, and the 
interactions between individuals are denoted by edges between nodes. For example, in a collaboration 
network, the scientists can be regarded as nodes and their joint papers can be treated as edges. From this 
collaboration network, we can find the collaborative relationship between any two scientists. However, many 
papers have more than two coauthors. Therefore, the collaborative relationship among any three or more 
scientists can not be described in the above simple collaboration network. Then, what kind of networks can 
solve this problem is an interesting issue, and the complex hyper-networks corresponding to hyper-graphs 
emerge in this context. 

In complex networks, an edge connects only two nodes, however, an hyper-edge of complex hyper- 
networks can contain more than two nodes, which can describe the multifaceted collaborative relationship 
appropriately. Though a hyper-edge can connect an arbitrary number of nodes, it is often useful to study 
hyper-networks where each hyper-edge connects the same number of nodes: a fc-uniform hyper-complex is 
a hyper-network in which each hyper-edge connects exactly k nodes. 

There has been a growing interest in the research of hyper-networks [10-14] recently. For example, in 
virtual enterprises, to respond to market changes rapidly and exploit unexpected business opportunities 
efficiently, the virtual breeding environment (VBE) actors collaborate, and share competencies, skills as well 
as resources [ll|. In order to describe the relationship among VBE actors, a hyper-graph (a hyper-network) 
is proposed as a meaningful logical structure, in which a hyper-path (i.e., a hyper-edge) represents the 



* Corresponding author. 

Email addresses: zhywu79810gmail.com (Zhaoyan Wu), jduanOipam.ucla.edu (Jinqiao Duan), xcfuOshu.edu.cn (Xinchu 

Fu) 



Preprint submitted to Elsevier 



November 14, 2012 



structure underlying a minimal cluster of enterprises. By the hyper-graph model, the formation process of 
a virtual enterprise can be presented and the pre-identified business opportunity may be caught. 

Recently, an evolving hyper-network model is introduced to describe real-life systems with the hyper- 
network characteristics [14] . Two evolving mechanisms with respect to a hyper-degree, namely hyper-edge 
growth and preferential attachment, are proposed to construct the hyper-network mode. The hyper-degree 
is defined as the number of the hyper-edge attached to that node. In this paper, a different evolving k- 
uniform hyper-network model is introduced, whose evolving mechanism involves the joint degree of k — 1 
nodes, which is defined as the number of the hyper-edges attached to the k — 1 nodes. 

On the other hand, synchronization as a typical collective dynamical behavior of complex networks has 
drawn considerable attentions recently [15-27]. However, less attention has been paid to synchronization 
of complex hyper- networks. Thus in this paper, synchronization of the 3-uniform hyper-network will be 
investigated for the first time. The rest of this paper is organized as follows. In Section 2, an evolving hyper- 
network model is introduced and several topological characteristics are studied. In Section 3, synchronization 
of a 3-uniform hyper-network is investigated and several synchronization criteria are obtained. In Section 4, 
several numerical examples are provided to verify the effectiveness of derived results. Conclusions are drawn 
in Section 5. 



2. An evolving hyper-network model 

An evolving mechanism with respect to joint degree to construct the /c-uniform hyper-network model is 
proposed in this section. In the scale- free network model put forward by Barabasi and Albert 0] , two simple 
evolving mechanisms, i.e., growth and preferential attachment, are proposed to construct the network model. 
Inspired by [2j, the following generation algorithm are proposed to construct the fc-uniform hyper- network 
model: 

(1) Growth: the hyper-network starting with mo (mo > k) nodes, at every time step, we add one node and 
m hyper-edges to the existing hyper-network, where m < C^" 1 = ( k _ 1 y ^_ k+1 y ■ 

(2) Preferential attachment: the probability IL; 1 i 2 ...i fc _ 1 of a new hyper-edge contains the new node and 
selected k — 1 nodes in the existing hyper-network depends on the joint degree of the k — 1 nodes at time t. 
Let di 1 i 2 ..,i k _ 1 denotes the joint degree of nodes i\, 12, ■ • ■ , ik-i, which is defined as the number of hyper-edges 
containing nodes i\, 12, ■ • • , ifc— 1, then 

-pr _ dj 1 i 2 ...j k _ 1 



d hh-ik-i 

l<il<i2<.--<ik-l<TO0+t-l 



where 1 < i\ < 12 < • • • < ik— 1 < mo + 1 — 1. 

Then, after t time steps, this constructs a fc-uniform hyper-network with mo+t nodes and rat hyper-edges. 
Define the hyper-degree d\ of node i as the number of the hyper-edges containing node i and node-node 
distance as the number of the hyper-edges in the shortest paths between two nodes. According to the 
generation algorithm and the relation between hyper-degree and joint degree, when a new node enter into 
the hyper- network, the probability of node i with hyper-degree d\ acquires a hyper-edge is 



n(4) = E 11- 

l<ii<---<j<---<ifc_i<mo+t — 1 

_ (fc-iK _ 

mkt 

Let i be the new node added to the hyper-network at time ti, and p(dh, U, t) be the probability that node i 
has hyper-degree d\ when it is being picked up at time t (t > U). Suppose that d l h is a continuous variable, 
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Figure 1: Hyper-degree distribution p(d/ l ) versus node hyper-degree on a logarithmic scale with the different parameters. 
The node number is N = 250. 



then probability U(d l h ) can be viewed as a continuous rate of change of d l h , i.e., 

where a is a constant and denotes the rate of change of d\. Since the new node i brings in m new hyper- 
edges, which has hyper-degree m at time tj, so the change of hyper-degree at tj is a — m. Solving equation 
([T]) with initial condition d l h (ti) = m, one has 

dl{t) = m( t -) k -^ 

H 

and 

, m . _k 

ti = (^)~t. 

d l 
u h 

Assuming that the time variables U have a uniform distribution with P(U) — 1/t, i.e., the new nodes are 
being added at equal time intervals, then one has 

P(d\(t)<d h ) =P{U>Q)^t) 

a h 

=i-p(u < (^)—t 
d h 

, m . _k 
=1 _ — ) — . 

d l 

Therefore, hyper-degree distribution P(d 1 h (t) — dh) can be derived as 

p{it) _ emm«u . ^a^a, 

which is the same as that of BA scale-free network when k — 2. Fig. [TJ shows the hyper-degree distribution 
of the 3-uniform hyper-network models with mo — 4, m — 1,3, 6 and t = 196, from which one can find 
that the 3-uniform hyper-network has power distribution property, i.e., scale-free property. Table I shows 
the average path length of the 3-uniform hyper- network models with mo = 4 and different pair of t and m, 
which exhibits small world property. 
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N = 50 


N = 100 


N = 150 


N = 200 


N = 250 


m 


= 1 


2.6176 


2.8893 


2.9459 


3.1752 


3.2785 


m 


= 2 


2.0539 


2.2786 


2.3688 


2.4823 


2.5124 


m 


= 4 


1.7683 


1.9461 


2.0135 


2.0849 


2.1437 



Table 1: The average path length of the 3- uniform hyper- network with mo = 4 and different pair of N and m, where N is the 
node number, i.e., N = t + itiq. 



3. Synchronization in hyper-network of coupled dynamical systems 

Complex network consisting of N identical linearly and diffusively coupled nodes can be described by 

N 

ii(t) = f( Xi (t))+e OijT(xj(t)-Xi(t)), i = l,2,--- ,N, (2) 

i=i 

where Xi(t) = (xn(t), Xi2(t), ■ ■ ■ ,Xi n (tj) T € R n is the state variable of node i, f : R n — !• R n represents the 
local dynamics of an isolated node, e > is coupling strength, V G R nxn is inner coupling matrix. Matrix 
A = (aij)jsrxN is outer coupling matrix, which denotes the network topology and is defined as, if there is a 
connection from node i to node j (i ^ j), then ^ 0, otherwise, ay = 0. According to the definition of 
the coupling matrix A, the network @ covers a variety of models ranging from simple weightless undirected 
networks to more complicated weighted directed networks. 

Similar to the network ([2"]>. a k- uniform hyper- network of N coupled systems can be described by 

N k 

it (*) = /(*<! (*)) + £ w-i^r^co-^Ct)), i 1 = l,2,.-. ,N, (3) 

*2,---,ifc = l j=2 

where i; (I = 1, 2, • • • , k) are different each other and a^^...^ is defined as: if there exists a hyper-edge 
contains nodes i\,i2,-- m , then ai 1 i 2 ...i h = 1, otherwise, a.^^...^. = 0. For fc = 3, the hyper-network ([3]) 
can be simplified as 

N N 

(t)+x k {t)-2xi(t)). (4) 

j=i k=i 

In the following, consider synchronization of the 3-uniform hyper-network which can be rewritten as 

N N N N 

Ht) = f(xi(t)) + eJ2(J2 ^k)T{x 3 {t) - Xi {t)) + a vk)r(x k {t) - x t {t)). (5) 

j = l k=l k=l j = l 

Let dij be the joint degree of nodes i and j (j ^ i), which is defined as the number of hyper-edges containing 
nodes i and j. Then one has 

N N N 

dij = Y a Vk = X! flifc J = " ' = Y afc J' i! 

fe=l fe=l fe=l 

and rewrite ([5]) as 

N N 

Xi(t) = f(xi(t)) + e dijT{xj(Jt) - Xi(t)) + e Y d ik T(x k (t) - Xi(t)). 
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Thus, 

N 

x i {t)=f(x i (t)) + 2s dijTixjW-Xiit)). (6) 

3 = ^d¥=i 



Define the diagonal elements of D as follows 

N N 



3 = 1,3^ 3 = ^,3^ 



Then, equation gives 

N 

Xi (t) = f(x i (t)) + 2eJ2dijTx j (t). (7) 
j=i 

In the subsequent studies, the hyper- network considered is always assumed to be connected, i.e., any pair of 
nodes is reachable along hyper-edge. It is easy to verify that the joint degree matrix D = (dij) is irreducible 
and its eigenvalues are = Ai>A2>--->Aat. Our objective here is to synchronize the network ((7|) with 
a given orbit s(i), i.e., 

xi(t) = X2{t) = ■ ■ ■ = xjyit) = s(t), as t — > oo, (8) 

where s(t) is a solution of an isolated node, namely, s(t) = f(s(t)). Here, s(t) can be an equilibrium point, 
a periodic orbit, or even a chaotic attractor. 

Let Xi(t) — s(t) + i]i(t) and linearize ([7]) about s(t). This leads to 

N 

fn(t) - Df(s(t))rH(t) + 2 £ ^rf y r % (t), (9) 

3=1 

where Df(s(t)) is is the Jacobian of / on s(t). Then refer to the proofs of lemmas 1 and 2 in (27j |. one has 
the following theorems. 

Theorem 1. Let = Ai > A2 > • • • > Xn be the eigenvalues of the joint degree matrix D — (dij). If the 
following N — 1 n- dimensional linear time-varying systems are exponentially stable 

i(t) = (Df(s(t)) + 2eA fc r)£W, k = 2, 3, • • • , N, (10) 

then the synchronized states |2P are exponentially stable. 

Proof. The proof will be given in the Appendix. □ 

Theorem 2. Suppose that there exists an n x n diagonal matrix P > and two constants S < and r > 
such that 

(Df(s(t)) + 5T) T P + P(Df(s(t)) + ~5T) < -rl n , (11) 
for all 5 < 6, where I n is an n x n identity matrix. If 

2e\ 2 < 6, (12) 

then the synchronized states are exponentially stable. 

Proof. The proof will be given in the Appendix. □ 
It is clear that the inequality (fT2j) is equivalent to 
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Figure 2: The second-largest eigenvalue of the joint degree matrix with m = 3. 

Therefore, synchronizability of the 3-uniform hyper-network (J5j with respect to a given coupling matrix can 
be characterized by the second-largest eigenvalue A2 of the joint degree matrix. 

Remark 1. Similar to the above discussions, synchronization of general fc-uniform hyper-network can be 
investigated as well. In fact, refer to the definition of joint degree, hyper-network ([3]) can be simplified as 

N 

Xi(t) = f(xi(t)) + (k- lje^dyr^t). 

4. Numerical simulations 

Firstly, consider a hyper-network consisting of 100 coupled Chua's oscillator (28[, which can be described 

by 

si \ / a(s 2 - si - /(si)) \ 

52 J = «1 - S 2 + S3 

53 J V -f3s 2 - 7S3 / 

where /(si) = bsi + 0.5(a — &) ( | si + 1| — |si — 1|) is a piecewise linear function, a > 0, /3 > 0, 7>0, a<0 
and 6 < 0. 

Fig. [2] shows the second-largest eigenvalue of the joint degree matrix D generating from the above 
evolving algorithm with mo = 4 and N = 100, where \2(N) is obtained by averaging the results of 10 runs. 
Then one has A2 < —4. In the numerical simulation, choose a = 10, j3 = 15, 7 = 0.0385, a = —1.27 and 
b = —0.68. Refer to the discussion in (2jJ, one can choose 8 — a such that the inequality (TTTj) holds. Then 
one can choose e = 0.5 such that the inequality ([TUl) holds. That is to say, the hyper-network generating 
from the evolving algorithm with m = 3 can achieve synchronization with any initial values. Fig. [3] shows 
the synchronization errors. 

Secondly, consider the synchronizability of 3-uniform hyper-network with different pair of m and N. Fig. 
U shows the second-largest eigenvalues of the joint degree matrix with different pair of N and m, where 
Aa(iV, m) is obtained by averaging the results of 10 runs. It is clear that \2(N, 1) > A 2 (^V, 2) > ••• > 
AziN, 6), which means the hyper-network with larger m has stronger synchronizability. Fig. [5] shows the 
synchronization errors with N = 100 and different coupling strength e = 0.3, e = 0.5 and e — 5, from which 
one can easily find that the hyper-network with m = 6 has stronger synchronizability. 

5. Conclusions 

An evolving mechanism with respect to joint degree is proposed to construct a fc-uniform hyper-network 
model. Based on a rate equation method, the hyper-degree distribution of the hyper-network is derived, 
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which obeys a power law distribution. Moreover, a complex hyper-network coupled with dynamical systems 
is introduced. By defining the joint degree of two nodes, the hyper-network can be simplified and its 
synchronization is investigated for the first time. Further, synchronizability of 3-uniform hyper-network 
with different pair of m and N is considered by calculating the second-largest eigenvalue A2 of the joint 
degree matrix. 
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Appendix 

Proof of Theorem 1: Let rj(t) = (r]i(t),T)2(t), ■ ■ ■ 7 i]N(t)) T G R Nxn , then equation © can be written as 



fl{t) = V (t)(Df(s(t))) T + 2eD V (t)T T . 
Choose an unitary matrix 4> = (cf>i, fa, ■ ■ ■ , (f>^) such that 

Mk = A fe( fe, k = l,2,--- ,N. 
Let A = diag(Ai, A 2 , • • • , \n) and 77(f) = $v(t) with v(t) = (vi(t),v 2 (t), ■ ■ ■ ,v N (t)) T G R Nxn , then one has 

v(t) = v(t){Df(s{t))f + 2eAv(t)T T 



Refer to the proof of Lemma 1 in Ref. [27], Ai = corresponds to the synchronization of the system states. 
Then if the N — 1 ?i-dimcnsional linear time- varying systems (|10D arc exponentially stable, the synchronized 
states © is exponentially stable. 

Proof of Theorem 2: Consider the following Lyapunov functions 



and 



v k (t) = (Df{s(t)) + 2e\ k T)v k (t), k = 1, 2, • 



,N. 



v k (t)=e(t)pm, k 



2,3,.-. ,N. 



Calculate the derivative of V(t) along (fTO)) gives 



Vk(t) = H T (t)((Df(s(t)) + 2e\ k T) T P 



+ 



P(Df(s(t)) + 2eX k T))at)- 



From inequalities (TTTT) and p^|) . one has 



Vk(jt) <e(t)((Df(s(t)) + 5TfP + P(Df(s(t)) + ST))t(t) 

<- remit) <~v k (t), 



where Ap is the largest eigenvalue of P, which implies 



Vkit) < ex P {--^}U(0), fc = 2,3,.-. ,N. 



S 



That is, the N — 1 time- varying systems (|10l) are exponentially stable, i.e., the synchronized states {5} is 
exponentially stable. 
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